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Abstract 

, We study the spectrum of fermionic modes on cosmic string loops. We find no fermionic 

, zero modes nor massive bound states — this imphes that vortons stabihzed by fermionic 

I currents do not exist. We have also studied kink-(anti)kink and vortex- (anti) vortex systems 



^3 



and find that all systems that have vanishing net topological charge do not support fermionic 
bound modes. 

1 Introduction 

Cosmic string loops have raised a lot of interest in the past decades. They are believed to form 

' in the evolution of cosmic string networks. While standard loops of cosmic string eventually 

decay, it has been suggested that fermionic or bosonic currents "living" on the loop of string can 

[ stabilize it against decay. This was first put forward in pj within a model of superconducting 

strings [2j. Such stable cosmic string loops are called vortons (see |3J and references therein). 

' In the case of bosonic currents, a scalar field condensate builds up in the core of the string and 

' leads to a conserved (Noether) current. Vortons have been constructed explicitly in a scalar 

O ' 

I field model with U{1) x U{1) symmetry [1]. 

' Fermionic currents appear due to the coupling of the string forming Higgs and/or gauge field 

. to fermions. The aim of this paper is the construction of the spectrum of both massless and 

^ ■ massive fermionic bound states on cosmic string loops. For an infinitely long straight string the 

spectrum is well known: the number of zero modes is given by an index theorem and equals 
the winding number of the string [Si |Bj , and there are massive bound states as well [71 [8] . The 
equivalent question for the string loop is hardly studied. An incomplete analysis of this case was 
done in [S1[TU]; they found that the number of zero modes is equal to the local winding number 
of the string. Although this seems to agree with the result as for the infinitely long straight 
string, we note that the straight string has a net winding whereas the loop does not. In [10] the 
existence of fermionic zero modes for a curved string was established analyzing small curvature 
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corrections to the straight string solution. However, it is not clear whether this analysis can 
be extended to string loops. In this paper, we asymptotically solve the fermionic equations of 
motion in the presence of a string loop. In contrast to [9], we will show that no fermionic zero 
modes or massive bounds states exists on cosmic string loops. And thus no fermionic vortons 
exists. 

In fact, we would not expect the existence of fermionic zero modes on string loops. The 
number of zero modes is an adiabatic invariant. One can adiabatically change the string loop 
to the trivial vacuum background, and the latter supports no bound states. Saying it in other 
words, the topological charge of the whole system is zero, and it follows from an index theorem [6] 
that there are no bound states. This statement can be verified explicitly in the much simpler case 
of a kink-antikink system. This system carries no net topological charge, it can be adiabatically 
deformed to the trivial vacuum. And indeed, solving the equations of motion, we find that 
the kink-antikink does not support fermionic bound states. The kink-antikink system is quite 
analogue to the set-up with two parallel straight strings with opposite winding number. For these 
simpler systems, we find that there are no fermionic zero modes for systems with vanishing net 
topological charge. 

Our paper is organized as follows. In section 2 we discuss the spectrum of fermion bound 
states on kinks and antikinks. In section 3 we move on to cosmic strings, and discuss single 
straight strings, two parallel strings as well as string loops. We end with conclusions in section 
4. 

2 Kinks and Antikinks 

We first examine the relatively simple system of kinks and antikinks. We expect that many of 
the results, which as we will see are topological in nature, carry over to the case of cosmic strings 
to be discussed in the next section. 

2.1 Bosonic background 

Consider a bosonic background of the form 



with r] > 0. B = corresponds to the well known kink solution of the (/(^-theory located at 
z = zi; the antikink solution is obtained by replacing r] —ij. B = —1 describes a kink-antikink 
configuration. This is only a solution of the field equations in the limit zi ^ oo. Finally, B = 1 
corresponds to a kink-kink solution. This set-up needs a bosonic potential with at least three 
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minima, and thus does not occur in the (/)^-theory. In principle, the bosonic potential can be 
reconstructed from the Bogomolny equations via dz4> = \/2V . 

The "thin wall" approximation, i.e. the limit that the domain wall is infinitely thin, is 
equivalent to the limit a — > oo. Eq. then reduces to 
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z < -zi, 

-Zi < Z < Zl, 

Z\ < z. 



(2) 



2.2 Zero mode solutions 

Consider a fermion with a Yukawa interaction £ 9 h^'i\)'\\). The Dirac equation in the background 
of kinks and antikinks described by ([T]) reads 



(3) 



To find the zero mode solutions we separate the longitudinal and transverse coordinates, and 
write 

■0 = a(t,2:,y)/3(z)e, (4) 

with ^ a constant spinor. Here a{t,x,y) solves the longitudinal Dirac equation and gives the 
dispersion relation. To find the zero mode solution we set E = 0, ki = for the moment (so 
that it is trivially solved), and a = 1; the case 7^ is discussed in the next subsection. To 
solve the transverse part of the Dirac equations we introduce the eigenspinors 
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The two projection eigenstates decouple in the Dirac equation, which becomes dzP± = ^h(p(3±. 
It has as solution 



P±{z) = N exp{^Bmjpz) cosh [a{z — zi)]^^"^*^'^ cosh [a{z + zi)] 



(6) 



with N a constant normalization constant, and = hrj the vacuum fermion mass. 

In the chiral basis for 7 matrices, see appendix IA.21 the eigenspinors of the projection 
operator ([5]) are of the form 
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There are four independent solutions = P±£,± with i = 1,2. They are normahzable if 



with 1pm labehng the 4 states V'i- 

Kink or antikink. To get the kink zero-mode solution we set i? = in dH]) to obtain 



The asymptotic behavior of the solutions is P± oc exp(=Fm^|2;|). The /3_ solution is non- 
normalizable, while /3+ corresponds to a zero mode solution localized at the kink at z = —z\. 
Two independent real solutions, or one complex solution, exist, constructed from '\\)\ = 
with i = 1,2. To find the solution for an antikink, we replace r] —rj. Now /?_ is the localized 
normahzable solution, while /3+ is non-normalizable. The profiles of the fermionic zero modes 
for kinks and antikinks respectively are shown in Fig. 1. 

Kink- kink. The kink-kink configuration corresponds to = 1. The asymptotic behav- 
ior of the solutions ^ then reads /?_ oc exp(m^|z|) which is non-normalizable, and /?+ oc 
exp(— 3m^|z|) which is normahzable. Both the bosonic background and the fermionic solutions 
are invariant under zi ^ —z\. This implies that we find only one complex zero mode solution, 
which is localized around z = —zi. The profile of the zero mode solution is shown in Fig. 1. 

Consider now the limit zi 0. The bosonic background then corresponds to a kink- 
configuration with (j) interpolating between —tj and 37]. The zero mode solution becomes /3± = 
exp(^m^2;) cosh[az]^^™''''/". As expected /?+ has the usual form of a normahzable zero mode 
on a kink The factor 2 difference in the exponent of the cosh is due to the fact that the 
topological charge Q oc (/)(oo) — (/>(— oo) is twice as big as for a usual kink. The factor exp{^m^x) 
comes from the fact that the B = 1 double kink solution is shifted horizontally by rj compared 
to the B = kink background, it is of the form cp = r] -\- 2r] tanh(a2;). 

Kink-antikink. The kink-antikink configuration corresponds to B = —1. Both zero mode so- 
lutions, although they are peaked at the kink or antikink, are non-normalizable since P±{zizoo) oc 
exp(m^|z|). Thus no normahzable zero mode solution exists. 

This result can be easily understood. The zero mode solution on a kink is a ^+-spinor, on an 
antikink a .^^-spinor. So pasting the kink and antikink together, one can never match the zero 
mode solutions at the origin because of the orthogonal spinors. This problem is independent of 
the distance between the kinks, in particular it does not disappear in the limit zi — > oo. And 
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P± = iVcosh [a{z + zi)]^™'*/'' . 



(9) 
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Figure 1: The profiles of the fermionic zero mode solutions /?+ (green) and (5- (blue) are shown 
together with the bosonic background (j) (red) for the kink (top left), antikink (top right), kink- 
kink (bottom left) and kink-antikink (bottom right) configurations, respectively. We have chosen 
zi = 10, 77 = /i = a = 1. 
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whereas the kink zero mode solution falls off exponentially at infinity oc e"^"^^ with (j) = rj, 
this no longer holds for the kink-antikink system. Indeed, extending the wave function past the 
antikink, cp flips sign, and /3+ cx e^^^^ blows up. 

Maybe in a cosmological setting the zero mode solution does not have to be extended all the 
way to infinity as there is a natural cutoff. The cutoff can for example be due to the presence of 
yet another kink or anti-kink (cf. the argument for the existence of a network of global strings), 
or (in a higher dimensional setting) the boundary of spacetime. We can compare the height 
of the wave function localized at the kink at z = —z\ with how it blows up at infinity if it is 
extended past the antikink at z = z\. This gives 

OO) ^va^{2z-z^)^ y ^) 

where we took the thin wall limit a oo. It follows that -Zi)| = \ip{z oo)| for z = 3zi. 
Thus the cutoff has to be smaller ^cutoff < 3^1 for the fermion to be localized on the kink at 
z = —z\. 

2.3 Massive bound states 

We now proceed to consider non-zero energy solutions. In the rest frame A: = 0, and the Dirac 
equation reads 

[i7°ao + i-i'd, -H\i) = 0. (11) 

Now ^!|_, which are eigenspinors of 7^ (see ([5]), ([7])), are not eigenspinors of 7'^. The fact that 
the energy is non-zero mixes the spinors i^±. 
We write the energy eigenstates as 

i,\ = ^-'''\fiX{z)t-iii'L{z)t) (12) 

with k = 1,2. The factor —i is put in for future convenience. and decouple in the 
Dirac equation. Their equations of motion are identical. In the following we will suppress the 
superscript, keeping in mind that the equations apply to both. The energy eigenstates satisfy 
-f^V'-B = -E'V'_B with 

H = idt = -P{-ifd,^H). (13) 
This gives the following set of coupled equations 

-E/3_ = {d, + H)(5+ (14) 

where we used 7'^Cit = 
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Let us analyze the asymptotic solutions. In the regions at spatial infinity, and for the case 
of multiple (anti)kink systems in between the (anti)kinks, the background field approaches a 
constant dz(j) = 0. In this limit the equations of motion decouple. Eq. (I14p can be written in 
the form 

dlP± = {m} - E^)p± (15) 



which has the solution 



P± = U±e^ f + D±e V / (16) 



with mf = h\(j)\ the vacuum fermion mass, and U±,D± integration constants. At spatial infinity 
one of the solutions is normalizable. Plugging back in (|14|) we find 



^ - ^ m?j — E'^ — hcj) J-) mj: — E"^ — h(j) 

TU " E ' d7 " E • ^^^^ 

In the core of a kink located at z = the bosonic background is lim^-^o </'(0) = V^lz + 0{z'^). 
To find the asymptotic solutions we write P± = h±z^^ near the origin, with n-t > to assure 
regularity at the core. The equations of motion (jl4p reduce to 

Eb+z'^+ = n_6_z"--^ -m/a6_z"-+\ 
-Eb.z''- = n+b+z''+^^ + mfab+z''++\ (18) 

The term proportional to nif is higher order. There is no general solution that solves the 
equations. There are however two special solutions. The first is the zero mode solution ([9]) with 
= 0, 6_ = 0, n+ = 0, which sets all lowest order terms to zero. The second special solution 
is a bound state solutions with E ^ and n_ = 0, n+ = 1. This choice allows the lowest order 
terms in the first equation in (llSp to cancel, whereas in the second equation the higher order 
terms come in at lowest order as well. 

Kink or antikink. The kink is known to have one bound state solution, which can be con- 
structed explicitly. For simplicity consider a kink located at zi = 0. To find the bound state so- 
lution we substitute /3+ = C(j)l3- in the Dirac equations (fUj) . Provided we choose C = —E/{ar]) 
and E'^ = arj{2h — a/rj), the equations are degenerate and read 

d-J^ = 7/tanh(az) - hj l3_ (19) 

which has as solution P-{z) = N cosh{az)^^^^"' with N a normalization constant. This bound 
state only exists for h > a/ {2rj). There is no bound state solution in the thin wall limit a/r] oo. 
Unlike the zero mode solution, the number of bound states is not an adiabatic invariant, it is 
not determined by an index theorem. 
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The above result agrees with the asymptotic analysis. Indeed, consider a kink located at 
z = —z\ in the thin wall approximation (see ([2]) with B = 0). In region I z S [—00 : 0] the bosonic 
background is = — and in region II 2; G [0 : 00] it is = r/. The normalizable solution in 
region i IS -i/)^ = e V ' ((7:j_^+ + (yi^_j, the solution m region ii is '(/'g = e v / + 
DU^-). Matching the solutions at 2; = requires 

UL/Ui = dL^/D^^ (20) 

which is impossible for E ^ 0. Setting = in the equations of motion, we see that the 
first equation of ()14p gives back the normalizable zero mode solution. Hence, in the thin wall 
approximation we find back the zero mode solution, but do not find additional massive bound 
states. 

Kink-kink Consider now the kink-kink system. Following the same strategy as for the kink 
system we write /?+ = C(pf3^. Now, however, we do not find a special solution to the equations 
of motion (jl4p . The difference with the kink case is the appearance of cross terms tanh(a(2; — 
zi)) tanh(a(z + zi)) which destroy the solution. One may argue that our ansatz for should 
be changed for the kink-kink system, but it is fixed by the requirement that at 2; = —zi the 
solution should approach the kink solution. 

In the limit — > the kink-kink system reduces to a single kink, which has a massive 
bound state. This bound state is very fragile. It dissapears in the thin wall limit. And from the 
kink-kink analysis it follows that it dissapears as well if the kink is deformed. 

The kink-kink system has one zero mode and no bound states in the spectrum. 

We have also tried to construct solutions to the equations numerically, but could not 
find solutions. 

Kink-antikink The story for the kink-antikink system is similar to the kink-kink system. We 
choose an Ansatz /?+ = C0/3_, which is dictated by the requirement that the solution approaches 
the kink solution at the kink at z = —Zc- No bound state solution is found, compared to the 
kink case it is the appearance of cross terms tanh(a(z — zi)) tanh(a(2; -|- zi)) that destroy the 
solution. 

The same conclusion also follows in the thin wall approximation ([2]) with B = —1. The 
solutions in the three distinct regions are given by The boundary conditions constrain 

the solution in the asymptotic regions: D± = in the z < —zi region, and U± = in the 
z > zi region. Matching the solutions at both z = —z\ and z = z\^ we find that the coefficients 
appearing in ()16p in the middle region —z\ < z < z\ satisfy Dj^/D^ = [/+/[/_ = 1, which is 
impossible to satisfy, even for E = Q. Hence, there are no bound states at all. 
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This can be understood as follows. The decaying solution at minus infinity is oc whereas 
the solution at plus infinity is oc and these orthogonal spinors cannot be matched in the 
middle. This spinor structure of the decaying solution is determined by the value (/)(iboo), hence 
is directly related to topological arguments. The decaying solution at plus and minus infinity 
is the same spinor iff the topological charge is non-zero: Q oc (j){oo) — (/){—oo) ^ 0. For the 
kink-antikink system, on the other hand (p{oo) — 0(— oo) = 0. The kink-antikink system can be 
adiabatically deformed to the trivial vacuum. As the vacuum carries no bound states, the kink 
anti-kink system has none either. 

The Dirac equation describes a fermion with a spatially varying mass term, i.e. a fermion 
living in an effective potential V^g = \h(j)\'^. For the kink-antikink system the effective potential is 
a double well potential. Experience with the analogue quantum mechanical double well systems 
suggests there should be two bound states, split in energy with the energy difference going to 
zero in the limit that the two wells are taken infinitely far apart. It seems from our results above 
that this intuition fails for the fermion zero modes. Why? 

The analogue with the quantum double well system breaks down. The double well potential 
in quantum mechanical problems can be adiabatically transformed by reducing the distance 
between the two wells, and for zi = one ends up with a single well. The number of bound 
states is preserved. What happens is that the two lowest bound states of the single well become 
the two energy split bound states of the double wells. Now taking zi ^ in the kink-antikink 
system, they will annihilate, and one ends up with the trivial vacuum. The vacuum does not 
support bound states, so it is not surprising that there are no zero modes on the kink-antikink 
system. The number of zero modes is given by an index theorem, it is {l/rj)((f){oo) — (j){—oo)), 
which is zero for a kink-antikink configuration. 

Note that our results disagree with those stated in [TTl [TU [T3] . In [T2] the energy is calculated 
for ^ = zizip-)B=o states, which is split in energy and goes to zero in the zi — > oo limit, just 
as in a double well potential. However, these are zero mode solutions to a background with only 
one kink or antikink, they do not satisfy the Dirac equation for the kink-antikink system. Using 
instead ^ = (-0+ it '0_)b=_i, the wave function is non-renormalizable and the energy blows up. 

3 Cosmic strings 

Having built up intuition on the existence of fermionic modes on relatively simple domain wall 
system, let us see whether the same kind of arguments apply to systems with cosmic strings. 

Consider two parallel strings, one with winding number (= topological charge) +n and one 
with winding number — n. There is no net winding around the string, and thus — by analogy 
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with domain walls — we do not expect fermionic zero modes in this system. The system can 
be deformed adiabatically to the trivial vacuum which does not support bound states. Looking 
at the spinor structure, zero modes on a string with n > (n < 0), have positive (negative) 
chirality. Pasting a string and "antistring" together, the zero modes cannot be matched in the 
middle because of the orthogonal spinor structure — in close analogy with the domain wall. 
This problem is not ameliorated in the limit that the strings are taken far apart. The solution 
that is localized at the string does not resolve the angular dependency at the antistring, and 
vice versa. It might be that a bound state with non-zero energy can cure this problem since it 
mixes the two chiralities. But we have seen that for the kink-antikink system this was not the 
case. 

The system of two parallel strings is closely related to the loop of cosmic string. In both 
cases there is no net winding number separating the system from the vacuum. Hence, we expect 
that the spectrum of fermionic bound states found on two anti-parallel strings (none!) is the 
same as for the string loop. This implies that there are no fermionic vortons, no loops of cosmic 
strings stabilized by fermionic currents. 

In this section we study fermionic bound states on straight strings and loops of string. In the 
next subsection, we discuss the bosonic string background. In section 13.2.11 and 13.31 we derive 
the fermionic spectrum for straight strings and a string loop respectively. To count the number 
of bound states an asymptotic analysis suffices. A numerical analysis is needed to obtain the 
profile of the bound states solutions, which can be done in the full system including back reaction 
effects. This involves systems of nonlinear ordinary and partial differential equations for straight 
strings and string loops, respectively; it is left as a future project. 

3.1 Bosonic string background 

In this section we discuss the asymptotic behavior of local cosmic string, both for straight strings 
and for loops. We neglect gravity effects, as well as back reaction effects from fermions (which 
will not effect the number of bound states). Moreover, we restrict ourselves to strings satisfying 
the Bogomolny limit. Note that in this limit there is no force between parallel straight strings. 
Details on the cylindrical and toroidal coordinates used can be found in Appendix El 
The flat space-time metric is 

ds2 = dt^ - hjdxl - hldxl - hldxl (21) 

The Lagrangian describing the gauged string (Nielsen-Olesen string) reads 

C = {D,cP){D^^r - \f,,F^^ _ ^ (|<^|2 _ ^2)2 ^ ^22) 
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with D^(f) = {d^ + ieA^)(j) and F^j^i, = d^Ai^ — diyA^. The static cosmic string configuration 
ahgned with the 3-axis, has a magnetic field confined to the string core and aUgned with the 
3-axis B = B3 = ^12/(^1/12)- The Higgs field winds around the X3 axis with winding number n 
(we take n> without loss of generality). The string energy per unit length is 



1^ ^ J V^~5^dxidx2 |(£>/i0)(D' 
= J ^/-9T<ixidx2 



2\2 



ho 



+ 



1 



F12 



2\2 



(23) 



where the integration is over the transverse coordinates xi, X2 which span the plane perpendicular 
to the string. The metric on this transverse plane satisfies ^—g^ = /ii^2- Following the 
arguments of Bogomolny — using partial integration and the identity [Di,D2](f) = ieFi2(f) — we 
write the tension as 




hi h2 



l\i\2 2\2 2 Fl2 

(101 -V) -er] 




F12 



h\ h 



1^2 



2 ■' ' hih2 

We have omitted the boundary term. The last term in the above expression is 

Q = -erf' J da;ida;2-Fi2 = -er/^ J ^ = "^T^nrf 



(24) 



(25) 



for a string with winding number n. Note that this is a topological number for a straight string, 
but not for a loop of string. For a loop, the string goes up and down through the transverse 
plane, and these contributions cancel. Hence, a string loop can decay, there is no topological 
barrier separating it from the vacuum. In the BPS limit e = A, the third term cancels, and it 
follows that n>Q. The minimum energy configuration with n = Q satisfies the BPS equations 



F12 
hih". 



hi h2 



-e(|</>p-r?2 



= 
= 



(26) 
(27) 



3.1.1 Straight string aligned with z-axis 



We use cylindrical coordinates with xi = r, X2 = ip and x-^ = z. The results for this case are 
well known, however, we repeat them here for completeness. The scale factors are hr = hz = 1 
and hin = r. Using the standard string Ansatz 



A, 



r?/(r)e-'^ 
n 

--a(r) 



(28) 
(29) 
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the BPS equations become 



t TL 

{dr + -D^)(t) = drf --{l-a)f = 

M^_e(|</,|2_r,2) = => n^-eV(l-/') = (30) 

with boundary conditions /(O) = a(0) = and /(oo) = a(oo) = 1. Note that the energy of the 
configurations remains finite, because the gauge field cancels the kinetic energy of the Higgs field 
(due to the phase dependence) at infinity. The magnetic field is given by B = V x A = drA^p. 
The solution is symmetric under rotations around the z-axis. 

First, we consider the equations in the limit r — > 0. The term proportional to a in the first 
equation tends to zero, and it follows / oc rl"L The second term in the second equation becomes 
r-independent and can cancel the first term for a oc r^, i.e. 




(forr^O) <; (31) 
a (X r 

Far away from the string core the fields tend to their vacuum values / ^ 1, a — > 1. To see how 
the fields approach these values we substitute / = 1 — cie~'^^^ and a = 1 — rc2e~'^^^ into the 
BPS equations. Upon taking the limit r — > oo we find C3 = V2er] = ruA the gauge boson mass 
(in the BPS limit uia = rrifj,) and 02/ ci = mA/n. We have: 

{f - 1 oc Q-'^AT 
(32) 
a — 1 oc (mAr-)e 



3.1.2 Cosmic string loop 

Consider a loop of cosmic string lying in the x-y-plane. We introduce toroidal coordinates with 
xi = V, X2 = u, = ^ such that the symmetry axis of the loop is the z-axis and the loop is 
aligned with the (^-direction. This set up is certainly less symmetric than that of the straight 
string since the fields now depend on the angle u winding around the string (in fact, to find 
the cosmic string loop, one has to solve partial differential equations in contrast to ordinary 
differential equations for the straight case). If the radius of the loop becomes very large, the 
u-dependence will practically disappear. The string Ansatz can be generalized to 

r?/(«,i;)e-« (33) 

ft 

--a{u,v) (34) 
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and we used the "radial" gauge Ay = 0. The boundary conditions are 

hm f,a = string core, 



D— »oo 



hm f,a = 1 spatial infinity, 

lim f,a = 1 origin. (35) 

(v,n)— >(0,7r) 

Note that we have assumed here that the loop radius is much larger than the string width 
K ^ iT^A^ J so that at the origin (and along the full z-axis) the bosonic fields approach their 
vacuum values. This requirement can be relaxed, all that is needed is that the fields are regular 
at the origin. However, in that case it is doubtful how well a cosmic string loop pictures the 
situation. 

The gauge field Au is non-zero to cancel the kinetic term of the Higgs field (j) at spatial 
infinity. The magnetic field is = —dyAu{u,v). The BPS equations are 

dyf + n(l - a)f = 0, duf = (36) 
nd^a + KKe^rf'il - f) = (37) 

We get two equations for / from splitting the real and imaginary parts. It follows that the Higgs 
profile function is u-independent; the only u dependence of the Higgs is in the winding phase. 
The Higgs profile is independent of the scale factors hi, and thus of the loop radius k. The gauge 
profile function a is generically u-dependent since h^hy depends on u. At infinity huhy is only 
a function of v and both / and a are functions of v only. Note that the u-independence of a far 
away from the string loop reflects the fact that there is no net winding around the string loop, 
no net topological charge. 

Let us consider the asymptotic solutions. For details of the behavior of the toroidal coordi- 
nates in these limits see Appendix [XI First consider the limit t> ^ oo in which the string core 
is approached and / ^ 0, a ^ 0. The BPS equations reduce to 

dyf + nf = 0, dya + const. • e'^^ = (38) 

which has as solution 



5r\'' _2y f Sr\ 



with 5r = \r — k,\ denoting the distance from the string core. The behavior is the same as in 
the core of a straight string (|3ip . The flelds approach their vacuum values far away from the 
string loop, corresponding to the limit u — > 0, u — > 0, and can be written in the form /=! — /, 
a = 1 — d/n. The BPS equations become 

dyf = d, dyd = huhym\f (40) 
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with niA = V^er]. At spatial infinity, corresponding to n = v — 0, the factor huhy tends to 
(k/v^)^. The equations are then solved to get (up to a normalization constant) 

/ - 1 (X ne-™-4^/^ ~ ne""*^'^ 

At the origin, corresponding to n — tt, u — >■ 0, the factor huh^ tends to (k/2)^. The solutions 
are up to a renormalization constant 

/ — 1 oc sinh ^ — — — ^ ~ sinh(m^r) ~ (rriyir) 

„ - 1 0, (^) (cosh - l) ^ (!!|^) (cosh(™..) - 1) ~ (m,rf (42) 

The approach to the vacuum is exponentially fast at spatial infinity (just as it is in the case of 
a straight string), but slower — power law — at the origin. The solutions are valid for w <C 1 
which corresponds to r 3> At at spatial infinity, and r <C At at the origin. The width of the string 



3.2 Fermionic spectrum 

Consider a Dirac spinor ^ = (^^), with 'ipL,'4'R l^ft and right handed two-component Weyl 
spinors. The fermion is charged under the U{1) of the string. The Lagrangian reads 

£ = 4a • DiPl + V-k • Di^R - ^H^yn ' V'rV'l (43) 

with D^ijji^R = (d^ + iqL.RC-A^)ipL^ii. Gauge invariance requires 1 — ql + Qr = 0, where we 
have normalized the Higgs charge to unity. For a Dirac spinor the two Weyl fcrmions ipL,i^R 
are independent, and can have different charges. For a Majorana fermion = —Qr = 1/2 and 
ijjR = icr'^tp'^ are charged conjugate of each other, so that the Majorana spinor satisfies the reality 
condition ipM = V'm- '^^^ reality condition decreases the degrees of freedom. The equations of 
motion are 

ia ■ Di/jl — X(l)ipR = 
iu-D^R-{\4>)*i,L = (44) 

For a Majorana the two equations are not independent, but complex conjugate of each other. 
3.2.1 Single straight string 

The transverse and longitudinal part of the solution decouple, and we can write 

il)L = a{z, t)PLir, (p)^L, tpR = a{z, t)pR{r, ip)^R. (45) 
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Here a and /Sl,/? are scalar functions of the transverse and longitudinal coordinates respectively. 
The constant Weyl spinors i,L,R are eigenvectors of the projection operator 



aV^a = i^, crV^eR = TCr, (46) 

with the upper sign for a vortex [n > 0) and lower sign for an anti- vortex (n < 0). Explicitly 
= (o) — i-i) n > 0, and = (5*) and = for n < 0, where the factor i 

in the right-handed spinors is for convenience. The Weyl spinors have the following properties: 
ia^(T^(7'^a^CL,R = -^L,R (chirality eigenstates), and a^'Ci = ie^'^^CR, (^""^R = -i&^'^^iL-, cr^^^L = 
^g±2(/3^^ and (j'^S.R = ^e^^'^ih- Using this the transverse part of the Dirac equation can be 
written as 

dr±i^^f3L = Ar?/(r)e±*(H-i)v'^^ 
{^dr^i^^dR = Ar?/(r)e^*(H-i)'.^^ (47) 

with /(r) the profile function of the Higgs field in (j29p . and as before the upper (lower) sign is 
for a vortex (anti- vortex). The phase-dependence can be resolved by choosing 

(3L,R{r, ^) = 6L,ii(r-)e*'^.«^, Il = Ir± (|n| - 1). (48) 

The equations reduce to 



(/l - nqiaY . 
<9r=F ]ol = ruffbR 



o j_{^R- nqRa)\ 

Or ± I bR = mffbL 



(49) 



where we used = r, and mj = \r] is the vacuum fermion mass. 

We can solve the Dirac equation asymptotically. For simplicity we restrict the analysis to 
the vortex solution with n > 0. At r ^ 00 the functions / — > 1, a ^ 1 and the 1/r term in the 
equations of motion is subdominant. The solution is 

hL,R = e^'^f' (50) 

There is one renormalizable solution. At the origin f ^ r"" and a ~ r^. We write bL,R ~ r^^^^. 
Then the equations of motion reduce to 

CRr"^-^ + lRr^R-^ = mfr^^+'' (51) 
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A singular solutions exist ior cl = II, cr = —Ir. Two regular solutions exist ior cl = II, 
cr = Il + u + I and cr = —Ir = n — 1 — II, cl = Il + 2?t.. They connect to the one normalizable 
solution at infinity. Both need to be normalizable at the origin 

j y^^d'^x^\ipL,R\'^ < oo j dr r\bL,R\'^ < oo 2cl,r + 1>-1. (52) 

This is satisfied for 

-l<lL<n. (53) 

Thus a string with winding number n has \n\ zero mode solutions with / = 0, .., \n\ — 1. The 
zero mode spinor is proportional to ^+ for a vortex with n > 0, and proportional to ^_ for an 
antivortex with n < 0. 

For |n| = 1 the zero-mode solution does not have any angular dependence. The Majorana 
solution (with ^pR = ifjf) can be given analytically [5]: 

The longitudinal part of the Dirac equation {a^do — (7^dz)ipL = becomes 

{do^dz)a = =^ a = a{t±z). (55) 

The zero mode on a string (antistring) moves at the speed of light in the minus (plus) direction 
along the string and has a dispersion relation E = \k\. 

3.2.2 Two parallel strings 

Consider now two parallel strings, one string located at x = xi and one string or antistring 
at X = —xi. If the distance between them is much larger than the string width, the bosonic 
background is well approximated by the two separated string solutions pasted together at x = 0. 

String-string. Consider two parallel strings which both have winding number n = +1. The 
bosonic background is then of the form (p = /ie*'^+/2e*'^ with /j — > at string i, and approaching 
the vacuum fi^l everywhere else. 

The zero mode solutions localized at the string at a; = —xi and the zero mode at the second 
string localized at x = xi have the same spinor structure. Moreover, near the string the angular 
dependence is resolved. Hence, they can be pasted together at x = 0. Indeed, we expect a 
zero-mode solution of the form ipR^L = (/?iC+ + P2^,+)r,l, with Pi localized at string i. The 
equations of motion read: 

dr + (/3i + /?2)l = Ar/(/i + hm + P2)r (56) 
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and a similar equation for the right handed spinors (interchanging L <^ R and a different sign in 
front of Dy,, see (I17|l ). The two solutions decouple and are approximately given by the solutions 
for the single string; since /i/32 — > at the location of the first string and /2/3i — > at the 
location of the second, the corrections are small. The important thing to note is that the zero 
mode solution localized at one string can be extended without problems to the region where the 
other string is, as the angular dependence is still resolved there — this is in sharp contrast with 
the string antistring system to be discussed shortly. 

This agrees with topological arguments. In the limit xi ^ the system reduces to a single 
string with winding number n = 2. Such a string has two zero mode solutions. As the strings 
are adiabatically separated, the number of zero modes remains constant. The two zero modes 
are the states localized at either one of the strings. This is different from the kink case. For 
strings, the number of fermionic zero modes equals the topological charge which is in Z, while 
for kinks the topological charge is in Z2 (asymptotic vacua either different or the same) and the 
number of zero modes does not exceed unity. 

String- Antistring. Consider now two parallel strings with opposite winding number, one 
string with n = +1 and the other with n = —1. The zero mode solution localized on the 
string and antistring have orthogonal spinor structure and decouple in the equations of motion. 
We can thus consider them separately. We write the zero- mode solution in the form 'ipL,R = 
(/3i^+ + f32^-)L,R, with Pi localized at string i. Let's concentrate on the Pi solution; the P2 
solutions are similar. The equation of motion for Pi is 



and a similar expression for L <^ R. The Pi solution does not solve the equations of motion 
near the antistring, where n = — 1 and the angular dependence is unresolved. Although /?i — > 
near the antistring, it is not equal to zero. Thus, the zero mode solution disappears from the 
spectrum. 

This agrees with the adiabatic argument. In the limit xi — > one ends up with the vacuum 
(n = 0). Since there is no topological charge, there are no zero modes. As the strings are slowly 
pulled apart the number of zero modes remains constant. 

The adiabatic argument should also hold for bound state solutions with E ^ 0. The non-zero 
energy mixes ^+ and ^_ spinors, which in principle might lead to a solution valid in the whole 
domain. We thus try the Ansatz 




(57) 




(58) 
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The equations of motion are then of the form: 



Ir 



m/(/ie''^ + /2e-*^)/32R 



m/(/ie''^ + /2e-*^)/3ifl 



(59) 



and a similar equation for L <^ R. To resolve the angular dependence we write /3j = e^^'^bi, with 
i = {1,2,3,4} = {1L,2L, 1R,2R}. Then the angular structure of the Dirac equation is of the 
form 



Set li = ^4 + 1 and I2 = h- This resolves the angular dependence of the first equation above. 
Inserting this into the second equation, one finds that no ^3 exists that can solve the angular 
dependence. Likewise, for li = I4 — 1 and /2 = ^4 + 2, the angular dependence for the first 
equation can be solved, but again no ^3 exists that solves the angular dependence of the second 
equation. There is no choice of li that resolves the angular dependence! And thus we conclude 
that the string antistring system does not support massive bound states either. 

3.3 Cosmic string loops 

In this section we discuss fermionic bound states on a loop of cosmic string. The total winding 
number of the string loop is zero, the bosonic configuration can be continuously deformed towards 
the trivial vacuum. And thus we do not expect any bound states to exist, in direct analogy with 
the kink antikink and vortex anti-vortex system. 

It is clear that the loop does not support zero mode solutions moving at the speed of light 
\14:\ I15j . Massless particles which go at the speed of light move along straight lines, and not 
along a curved path. (Indeed, along a curved path k ^ and a dispersion relation E"^ = k'^ is not 
possible). But a priori this does not exclude a zero mode solution in the limit /c ^ 0. In Ref. [9j 
the existence of zero mode solutions was claimed. However, their analysis is incomplete. They 
failed to check the behavior of the zero mode solutions at the center of the string, where the 
solutions is singular. We show this explicitly in the appendix [Bl where we give the asymptotic 
solutions of the Dirac equation. 

The obstruction to bound states is the same as for the vortex anti-vortex system. For a loop 
lying in the (a;, y)-plane, consider the solution in, say, the x = plane. The angular dependence 

^The single straight string is known to have massive bound states [3 [8]. These states do resolve the angular 
dependence in the Dirac equations. 



_|_ g«('2-l)¥' 



gi(/4 + l)<^> _j_ gj(/4-l)</3 
gi(/3+l)i/p _|_ ^i{h-l)<^ 



(60) 
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of the equation of motion is exactly the same as for the vortex-antivortex. Since in the latter 
case the angular dependence could not be resolve, the same is true for the string loop. There 
are no bound state solution. 

To see this all more explicitly, let's look at the Dirac equation. To describe bound states on 
loops of cosmic string we try the same strategy as before and split the solution in a longitudinal 
and transverse part 

V'L,/? = a{t, ip)(3L,R{u, v)Cl,r, (61) 
with S^L^ji eigenstates of the projection operator 

aV'^a = Hl a'^a'^^n = (62) 
Explicitly, the eigenspinors corresponding to the upper sign are 

and the eigenspinors corresponding to the lower sign are 

In contrast with the straight string, the eigenspinors are not constant but depend on the angle 
if along the string loop. The ip dependence is needed to resolve the (/j-dependence in the Dirac 
equation. For future use, we note that the spinors satisfy a^^L,R = —£,l,r, d^p^L = '^\^R) and 
d^iR = T^ih- 

The longitudinal part of the Dirac equation is [dt ~ '^^l^ = and a similar equation 
for ijjR. Now since ^^pE,± / an extra term appears and the equation reads 

dta T a + ^a^R = (65) 

r J 2r 

This cannot be solved as for the straight string. First of all, the equation depends through 
= r{u, v) on the longitudinal coordinates, and the separation of variables is no longer valid. 

Secondly, there is an extra term since the eigenvectors are 99-dependent, but one could hope that 

this term can be "absorbed" in the longitudinal part of the Dirac equation. 

Nevertheless, let us try to split the frequency part and see where we get. Try the Ansatz 

'^L,R = e~*'^*(/3i^+ + Pi^-)r^l, with {C±)l,r the spinor eigenstates of ([B^ . The transverse Dirac 

equation in cylindrical coordinates is 

iEl3iL + {Dr-iD, + ^)l32L + iX(^l32R = 

iE^2L + {Dr + iD, + ^)PiL-iX<P^iR = (66) 
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and a similar equation for L ^ R. Now transform to toroidal coordinates, see appendix [Xl for 
details. Using Dj- + iD^ = —{2i/a) sin^ x(-^u + iD^) and l/(2r) = sin x sinx/ (asinh v) we get 

( 9^/^ isin?sin£\ ^ a^(t> ^ 

^?ii;/3..+ fsm^eP. + .Z?.)-^|^>i. + ^/?i« = (67) 
2 \ isinnw J 2 

If we follow the same strategy as for the straight string, we now want to remove the u depen- 
dence, i.e. the dependence on the coordinate winding around the string. This, however, seems 
impossible. Remember that X) X depend on u, so sinx contains a e™/^ and a e"*"/^ term. The 
phase structure of the equation above thus is 

PiL + (e'" + 1 + e-^")/32L = e^""/?2R (68) 

and similar for 1 <^ 2. There is no way to remove the phase dependence. The situation is 
completely analogous to the vortex antivortex system, whether E is non-zero or not. 

4 Conclusions 

In this paper we have studied the fermionic spectrum on different topological defects including 
systems of soliton-antisoliton systems. We find that the number of zero modes is given by 
topological arguments, just as in the case of a single soliton. Any configuration that can be 
continuously deformed to the trivial vacuum does not support fermionic bound states. Although 
this sounds like a trivial conclusion, statements to the contrary can be found in the literature. 

To built up intuition, we first studied the simpler kink-antikink systems. While a kink, 
antikink and kink-kink system all support one fermionic zero mode, a kink-antikink has none 
due to the different spinor structures of the fermionic mode associated with the kink and antikink. 
This is in complete agreement with topological arguments. Moreover, we find that no massive 
fermionic states on any system of kinks and anti-kinks exists. 

Cosmic strings with winding number n are known to support n fermionic zero modes. Ac- 
cordingly, we find that a system of two parallel BPS strings with n = 1 each supports two 
fermionic zero modes, while a string-antistring pair has no zero modes. In addition we find that 
a loop of cosmic string (which can be approximated by a string-antistring pair) has no fermionic 
zero modes. Note that this contrasts with the results in [9} llOj. 

Cosmic string loops stabilized by fermionic or bosonic currents, so-called vortons, are believed 
to be important for cosmology. Our results indicate that cosmic string loops with fermionic 
currents on them do not exist. 
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Of course, we have neglected the back reaction of the fermions on the bosonic fields and have 
only studied large loops. We believe, however, that taking these two points into account will 
not change our main conclusions. 
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A Cylindrical and toroidal coordinates 



We use the flat space-time metric with signature (+,—,—,—). Consider the three spatial di- 
mensions with a metric of the form (1211) 



ds^ = -hjdxj - hldxl - hldxj. (69) 

In Cartesian coordinates {xi,X2,X3) = {x,y,z) with hi = 1. Cylindrical coordinates are most 
convenient to describe a straight string aligned with the z-axis. Then (xi, X2, xs) = (r, ip, z) and 
hr = hz = l,h^ = r. The relation to Cartesian coordinates is {x,y,z) = {r cos (p, r simp, z) . The 
Pauli matrices in cylindrical coordinates are: 

Toroidal coordinates are most convenient to describe a loop of string lying in the x-y-plane. 
Then (xi,X2,X3) = {v,u,ip) with 

h — h - 



cosh V — COS u 
K sinh V 

K = Z • (71 

cosh V — cos u 



The relation with Cartesian coordinates is 



, K sinh u cos 09 k sinh u sin w Ksinn 

{x,y,z) = {^ ^, — ^ . 72 

cosh V — cos u cosh v — cos u cosh v — cos u 

At the core of the string the radius of the loop is k, and its arc length ndip with ip G [0, 27r). The 
coordinate u E [0, 2tt) winds around the string core, whereas v gives the radial distance away 
from the core (but note, for f — > oo the string core is approached, whereas f — > corresponds 
to far away from the core). The Pauli matrices in toroidal coordinates read: 

^ 11 sinusinht; e~*'^(cosncosht; — 1)\ 

<7 = r . ) (73) 

cosh V — cos u ye*"^ (cos n cosh t> — 1) — sin tt sinh v / 

1 / — cos ti cosh t> + 1 e"**^ sin n sinh , ^ 

<y = r . , (74) 

cosn V — cos u \ qI-v gjj^ ^ ginh v cos u cosh v — II 

= [ . . 75 

J 

We also introduce the notation 

X = {u + w)/2. (76) 
22 



The equations of motion may be simplified using sin X sin % = i (cosh v — cos u) and sin^ x = 
^ (1 — cos u cosh V — isinu sinh v) . 
For 1-forms we use the notation 

A = Y,Adx' = ^e' (77) 
tli 

with vierbein = /ijdx*. For example, for a straight string in cylindrical coordinates 



r 

^^||dx + ^v^dy = -A^^dx + A^^dy = A^dx + Aydy. (78) 



And similarly for higher forms: F = Fi2dx^dx^ = /^^e^e^. 
A.l Asymptotics 

The core of the straight string aligned with the z-axis is approached in the limit r = \Jx^ -\-y^ 
(and ill defined), whereas asymptotically far corresponds to r — oo. 

In toroidal coordinates the core of the string is approached in the limit w — ^ oo (and u ill 
defined). Then {x,y,z) — (kcos Ksinc^, Ksinue"^) and {hy,hu,hip) — {Ke'~'" , Ke~'" , k) . The 
core of the string loop is at r = \J x^ -Vy^ = k- Write r = k + 5r, then we can express the limit 
V — ^ oo in terms of (5r — ^ 0. More precisely 



V = — arccosh 

and thus 



lim e-" ~ hm — (80) 
(r, K have dimension of length, while u, v, (p are dimensionless.) 

Spatial infinity corresponds to u = v ^ 0. In this limit coshw — costi ~ f^, sinhf v and 
thus 

{x,y,z) — > {-cos(p,-sm(p,-) 

V V V 

{K,hu,h^) ( _) (81) 

and 

lim V oc lim — oc lim — (82) 

with R = \J x^ -\- y^ -\- = Vr^ + oc r. In the limit k <^ R, the distance from the string core 
approaches the distance from the origin R. 
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Finally, the limit u = tt,v ^ corresponds to the origin, the center of the loop. Indeed, in 
this limit 



and 



{x,y,z) — > {-V cos If, -vsmip,0} 
iK,hu,h^) (83) 



lim V oc lim — . (84) 

It— ►7r,ti"»0 r—rO K 



A. 2 Notation and convention 

Most of the solutions are given in terms of eigenspinors, and are independent of the basis of 
7-matrices. When we give explicit solutions, we use the chiral basis: 

with o"'^ = (1,CT*), (T^ = (1, — (T*), and cj* the Pauli matrices. 



B Asymptotic solutions of [j9| 

In this appendix we show explicitly that there are no fermionic zero modes on a circular loop of 
cosmic string. This is in contrast with the results in [9] . Following [9] we remove half of degrees 
of freedom by considering Majorana fermions with ippi = V'l ~ ^'^'^l- The equation to solve is 
then (|67p with E = 0. Using the Ansatz 

f3 = Xe^'(^+*) + y*e^("-')(^+^) (86) 

this gives the two coupled differential equations 

( ■ -i - ( r. , ^ \\ , «sinxsinx\ „ ^mf 

sm Y [Oy + a) H X = Y 

V V ^ ^ 2 V 2sinhw J 2 

with m = Xt] the vacuum fermion mass. These equations agree with those found in [9]. 



Infinity. Let us first consider the behavior at infinity and write 



X = l{u + iv) = -^e''^. (88) 
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Spatial infinity is approached in the hmit u = u — > 0, or equivalently t ^ oo,<p ^ 7r/4. Further, 
the gauge field and Higgs field approach their vacuum values and thus a{v) 1 and f{v) 1. 
Then using the approximation sin^ ~ X ^iid after redefining the fields 



871) becomes 



X = aJ^, Y = bJ^, (89) 

y x-x y x-x 



Km 

XX [Ox -ip)A = —B 

XX{dx-w)B = '^A (90) 
where we defined p = n/2 — I. Combining the two equations gives 

— j A. (91) 
Using (I88|) we write this in terms of t, ip to take the limit r — > oo. The partial derivatives become 

dx = —(dt + f^), d^ = -^{dt--d^) (92) 

and XX = t^"^ ■ Using these expressions in (|9T]l and keeping only the terms with highest power 
in t gives: 



\d^,A - ^dtd^A - f A + 0(t-2) = (93) 



A ' At ' V 2 

from which it follows that 

hm A = Ce-"""^ oc e^""^ (94) 

r— >oo 

with C a normalization constant, and where we used t cx oc n/r. There is a similar solution 
for B = &~™I'^A. Putting back all the factors we get for the normalizable solution: 

= ^^^lV2/t ^ g^./4g-^(n~0^/2/^) . (95) 

Loop core. The core of the loop is approached in the limit f — > oo or x — > iv/2 — > zoo. In 
this limit sinx — > ie^/^. Further /, a — > at the core. Taking this limit in the Weyl equation 
(|87|) ■ only the first term survives and gives 

{d^ + il)X = Q X~/x(x)e-''^ 
(d^-{n-l))Y = r ~ /y(x)e*("-')^ (96) 

Then 

f3 = /x(x)e*'^ + /y(x)e*("-')^ 

Iv {■n, — l)v 

= fx{v)e~^ + fY{v)e 2 for u ^ oo (97) 

Since the functions fx{v), fviv) are unconstrained, nothing much can be said about the number 
of solutions at the string core. 
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Origin. Finally consider the zero mode solution at the origin (u, v) (it, 0) or (x, x) ~^ 
(7r/2,7r/2), where / ^ 1, a ^ 1. Note that this case was not studied in 0. The Weyl equation 



871) becomes 



idx-ip) + lT7—^]y = —X (98) 



1 \y '^"^ 

where we have used the abbreviation p = ^ — I. 

The problem is already apparent here: X ~ X ~^ the second term on the left hand side 
of both equations blows up. Redefining the fields 



the equations simplify to 



X = A. / Y = B. / (99) 

\ x-x y x-x 



{d^-ip)B = (100) 



2 

which can be combined into a single equation for A 



, , 2 

2 " 



(^0^_,p)^d._ip)A=(^) A (101) 



and likewise for B. Explicitly, we have: 

d^d^A - ip{d^ + d^)A - [p^ + {nm/lf) A = {). (102) 

We will try two type of Ansatze. First, we choose A = A[x + x) = ^(^)) which implies 
A" = 2ipA' + {p^ + (Km/2)2)A and the prime denotes the derivative with respect to u. This has 
as solution 

A^C^{±nm/2+^p)u «-_^-0 ^^^^^^ ^^^^^ 

Further B (x A where the constant of proportionality is given by the equations. 

The second type of Ansatz is A = A{—i(x ~ x)) = ^(^)) which gives A" = {p^ + {K'm/2)'^ A 
where the prime denotes derivative with respect to v. This has as solution 

A = Ce±Vp^+(W2)2. u^^-.o ^^^^^ ^^Q^^ 

Both solutions do not tend to zero at the origin. On the contrary, plugging the above expressions 
back in, we find that zero mode solution 

f3=^J^ (^e*'" + i?*e'("-')'') "^^"^V l/v (105) 
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Although the solution blows up in the origin as 1/v the integral J dvduy/g2il)'^ip is finite. Nev- 
ertheless, this is not a valid solution. The spinors (I63p . ()64p are only well defined at the origin 
if /? 0. Note that a similar argument holds for the full z-axis with t; = and n 7^ 0. The 
solution is thus singular on the symmetry- axis. Non-zero energy will affect the exponents in 
B, but not the dominant behavior of f3 in the limit u ^ vr, t; — > 0. And thus non-zero energy 
cannot alleviate this problem. 
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